Photoinduced time-resolved electrodynamics of superconducting metals and alloys 
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The photoexcited state in superconducting metals and alloys was studied via pump-probe spec- 
troscopy. A pulsed Ti:sapphire laser was used to create the non-equilibrium state and the far-infrared 
pulses of a synchrotron storage ring, to which the laser is synchronized, measured the changes in 
the material optical properties. Both the time- and frequency- dependent photoinduced spectra 
of Pb, Nb, NbN, Nbo.5Tio.5N, and Pbo. 75610.25 superconducting thin films were measured in the 
low-fluence regime. The time dependent data establish the regions where the relaxation rate is dom- 
inated either by the phonon escape time (phonon bottleneck effect) or by the intrinsic quasiparticle 
recombination time. The photoinduced spectra measure directly the reduction of the superconduct- 
ing gap due to an excess number of quasiparticles created by the short laser pulses. This gap shift 
allows us to establish the temperature range over which the low fluence approximation is valid. 

PACS numbers: 74.40.-|-k, 74.70.Ad, 74.25.Gz, 78.47.-l-p 



I. INTRODUCTION 

When a superconductor is subjected to a pulse of 
above-bandgap electromagnetic radiation, excess quasi- 
particles (QPs) are created. The number of excess QPs 
created, and the way in which they recombine back into 
Cooper pairs, are processes of fundamental importance. 
Burstein and co-workers- made the first attempt to es- 
timate this recombination rate. Assuming a radiative 
process (emission of a photon) they calculated a recom- 
bination time around 0.4 s. Ginsberg's^ tunneUng data 
in a Pb film, into which excess QPs had been injected, 
showed recombination times many orders of magnitude 
smaller (10"'' s) giving indications that the process was 
dominated by emission of phonons rather than photons. 
Subsequent calculations^^i^i^ confirmed the faster re- 
combination rate when phonon emission is considered. 
However, Rothwarf and Taylor^ showed that the mea- 
sured (or effective) lifetime was longer than the "bare" 
or intrinsic recombination time due to phonon-trapping 
effects. 

Many groups have studied the problem of nonequilib- 
rium superconductivity by injection of unpaired quasi- 



particles in tunnel j unctions i^iSiSiifliii In these steady- 
state measurements, the excited state is detected through 
an average response that deviates from equiUbrium and 
the effective pair relaxation time is inferred, assuming 
that the process has a simple exponential decay. 

As an alternative, the non-equilibrium state can be 
generated by breaking Cooper pairs. Testardii^ showed 
that light of enough intensity could destroy the supercon- 
ductivity in a process not attributable to lattice heating. 
In this process, photons with energy larger than 2A (A 
being the superconducting gap) break pairs. Parker and 
Williamsi2i used this phenomenon to measure pair recom- 
bination rates in illuminated tunnel junctions. 

Owen and Scalapinc^'* described the non-equilibrium 
state by an effective chemical potential (/j,*) and showed 
that the excited state is characterized by a shift in the 
superconducting gap. Sai-Halasz et al.^^ used a pulsed 
laser to create the excess number of QPs and analyzed 
the changes in the microwave reflectivity. Their data sup- 
port the theory of Owen and Scalapino for the low-flucnce 
regime, where the number of QPs created by the photon 
absorption is much smaller than the number of unpaired 
QPs produced by finite sample temperature. The gap 



2 



shift in the photoexcited state was also detected by Hu 
et a/.A§. in tunnel junctions; this experiment, which had 
a time resolution of 30 ns, was the first direct measure- 
ment of the relaxation process. Parker then proposed 
that the photoexcited state could be described by BCS 
theory with an effective temperature (T*) based on the 
phonons with energy greater than 2A,^'' For small per- 
turbations both T* and fi* models are indistinguishable. 

The major drawback of most previous work is the indi- 
rect manner by which the relaxation rates were measured. 
JohnsonliSi was the first to use fast lasers and electronics 
to look at voltage transients with time resolution in the 
range of 100 ps. He directly measured the effective relax- 
ation times in Nb films and found a relaxation composed 
of a double exponential decay, suggested to be intrinsic. 
Federici et al}^ used coherent THz techniques to look at 
the pair- breaking dynamics, finding that pairs are broken 
in picoseconds. 

More recently, all optical ultra-fast pump-probe tech- 
niques were used to study the non-equilibrium dynamics 
of high-Tc superconductorsi^°i^^i^^i^^^'^i^^i^^i^'^ In these 
materials, relaxation times were three orders of mag- 
nitude faster than in metallic superconductors, but no 
direct gap shift was measured. Although the gap sym- 
metry of the high- Tc materials plays a major role in the 
relaxation times, qualitatively the process is similar to 
the one observed in BCS superconductors.^^ All optical 
pump-probe spectroscopy also proved to be very useful 
in the understanding of the non-equilibrium dynamics of 
the double gapped MgB2.^2*^ii2S 

In this paper we extend our far-infrared pump-probe 
spectroscopy work done on Pb thin filmsH to Nb, NbN, 
Nbo.5Tio.5N and Pbo.75Bio.25 and we develop several is- 
sues linked to the recombination process. This paper 
is divided as follows: in Sec. ^ we describe our samples, 
the experimental setup and the characterization methods 
used. In Sec. IIIII we discuss the equilibrium far-infrared 
spectra of our samples. Section IIVI deals with the non- 
equilibrium excess QP state while our conclusions are 
presented in Sec. [3 In the Appendix we discuss the con- 
sequences of sample heating by the laser and methods for 
minimizing this effect. 



II. EXPERIMENTAL 

A. Samples 

Our samples were thin metallic films of Pb, 
Pbo.75Bio.25, Nb, NbN, and Nbo.5Tio.5N. AU these metals 
were deposited on 0.5 mm thick c-axis sapphire, except 
for the NbN, which was deposited on 0.5 mm thick (100) 
MgO crystal. The film thickness were not directly deter- 
mined. As it will be discussed in Sec. IIIII the important 
quantity is the sheet resistance i?n which can be calcu- 
lated from the normal state transmission. From i?Q we 
can estimate roughly the thicknesses of our films, obtain- 
ing values ranging from 50 to 200 A. Several films of each 



composition were grown. 

Generally, the far-infrared transmission decreases ex- 
ponentially with the film thickness, therefore decreasing 
the signal to noise ratio. However, Tc is usually higher for 
thicker films. In this paper we show in detail data for the 
films with the highest T^s consistent with a normal-state 
far-infrared transmission between 10% and 30%. 

Pb and Pbo. 25619.75 films were evaporated in situ at 
80 K and the samples were kept below 100 K throughout 
the measurements to avoid annealing and normal state 
conductivity change. The transition temperatures were 
optically determined (see 'Methods' below) by the in- 
crease in the far-infrared transmission that begins at Tc- 
Twenty Pb films were deposited and we measured Tc be- 
tween 7 and 7.2 K. The four Pbo.25Bio.75 films showed the 
superconducting transition at 8.0 K. The Tc's in these 
samples were independent of the deposition time and 
hence of the film thickness. 

Nb films were sputtered and kept in vacuum to avoid 
oxidation and degradation of the superconducting prop- 
erties. Two different samples showed (optically deter- 
mined) Tc of 5.9 K and 5.4 K. We show here the data for 
the latter. NbN was produced by reactive sputtering of 
Nb in N2 atmosphere.2i2^ The sample had Tc = 13.5 K. 
Nbo.5Tio.5N films were sputtered in an N2 atmosphere. 
For thick films (in the /im range) this technique produces 
samples with Tc ~ 14.5 K. Our films were much thinner 
(100-200 A) and showed suppressed Tc's, below 11 K. 



B. Apparatus 

In our measurements we utilized the far-infrared pho- 
tons from beamline U12IR in the VUV synchrotron stor- 
age ring of the National Synchrotron Light Source at 
BNL4& 

Spectrally-resolved measurements were recorded from 
5 to 60 cm~^ with a lamellar grating interferometer'^^ for 
the Pb, Pbo.25Bio.75 and Nbo.5Tio.5N films. A Bruker 
IFS 66v was utilized to measure Nb and NbN samples 
from 5 to 100 cm"-'^. Most of the data were obtained with 
a 1.3 K bolometer detector although a few measurements 
were taken with a 4.2 K bolometer. Black polyethy- 
lene and fluorogold filters gave an upper frequency limit 
of 100 cm~^. The equilibrium transmittance spectra of 
NbN was measured at ESPCI in a Bruker IFS 66v using 
a Hg-arc source and a home-made cryostat. 

Low-temperature transmittance data for the Pb and 
Pbo.25Bio.75 films were measured using a ARS He-Tran 
cold finger cryostat. The samples (10 x 10 mni'^) were 
solidly clamped to the cold finger using indium gaskets. 
A 5 mm diameter circular aperture was used, and the 
remainder of each sample remained in contact with the 
cold surface. The NbN and Nbo.5TiQ.5N films were mea- 
sured in an Oxford Optistat cryostat, with the samples 
immersed in cold gas or liquid helium. The Nb films 
were measured in both cryostats. In some experiments, 
we measured the temperature-dependence of the overall 
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far-infrared transmission (without frequency discrimina- 
tion by a spectrometer) through the sample. In this case 
we modulated the beam at around 100 Hz with an optical 
chopper and used a lock-in amplifier. 

To create the non-equilibrium photoinduced state, a 
mode-locked Ti:Sapphire laser was synchronized to the 
rf cavity of the VUV synchrotron storage ring. The laser 
produces a maximum of 1 W average power distributed 
in pulses with a repetition rate of 52 MHz, the same repe- 
tition rate as the synchrotron pulses. The laser pulses are 
approximatelly 2-3 ps long at a photon energy of 1.5 eV. 
The laser is used to produce excess quasiparticles in the 
superconductor and the far-infrared synchrotron pulses 
probe the changes in the sample transmission caused by 
these broken pairs. The time resolution is determined by 
the synchrotron pulse width, and is ~ 200 ps. The inter- 
pulse separation limits the window for relaxation to 19 
ns. More details on the facility can be found in Ref. IsM 
Typical laser powers used in our experiments were in the 
range of 10 to 100 mW (0.2 to 2 nJ per pulse) evenly 
distributed over the 5 mm clear aperture on the samples. 
Such low power is required to minimize heating and to 
keep the material in the low-fluence regime. 



C. Methods 

Two types of transmission were employed: (i) 
we recorded the overall far-infrared light transmitted 
through the sample without energy discrimination and 
(ii) we measured frequency-resolved spectra between ~ 5 
and 100 cm~^(60 cm~^ for measurements in the lamellar 
interferometer) . The overall transmitted light is obtained 
simply by placing the sample in the optical beam between 
the source and the detector. This overall transmission 
is the integrated sample transmission below 100 cm^^, 
weighted by the emission spectrum of the synchrotron 
and by the transmission of the filters, the cryostat win- 
dows, and the low-pass filter in the bolometer cryostat. 
So long as the cut-off frequency is higher than about 3A 
(which is our case), the transmission of a superconduc- 
tor increases in the superconducting state. This prop- 
erty can be used to measure the sample Tc- Figure Qfa) 
compares the increase in the overall far-infrared transmis- 
sion to the magnetization measured in a Quantum Design 
MPMS squid magnetometer for one of our samples. The 
optically-determined is about 0.5 K higher than the 
value obtained from the magnetic measurement. This 
discrepancy is attributed to the fact that the far-infrared 
transmission senses resistive transitions, and percolation 
of one superconducting path may happen before the en- 
tire sample becomes superconducting. As shown in Fig. 
^b) this overall transmission can also be easily related 
to the number of Cooper pairs in the system. The curve 
in Fig. ^b) is not universal and must be determined for 
each sample measured. Note that it is approximately but 
not exactly linear. 

The non-equilibrium state created by the laser pulses 
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FIG. 1: (Color online) (a) The left axis shows the increase 
in the far-infrared transmission due to the superconducting 
transition in one of our Nbo.5Tio.5N samples. The solid line 
is a BCS calculation for the temperature evolution of the pen- 
etration depth A'^(0)/A^(r), a quantity which is proportional 
to the number of pairs in the system. On the right axis we 
plot the magnetization of the same sample. Note the sup- 
pressed zero, (b) Overall far-infrared transmission vs. calcu- 
lated A^(0)/A^(r). (c) Differential transmission measurement 
(o) and its integral (□) for a Nbo.5Tio.5N film at 2 K. Point 
'A' in the differential curve corresponds to the maximum up 
slope in the integrated curve and point 'B' to its maximum 
down slope. 



was, for the most part, studied using the overall trans- 
mitted light through the sample. These data allowed us 
to measure the variation of the transmission as a function 
of the pump and probe delay, yielding the temperature 
dependence of the amplitude of the photoinduced signal 
and its relaxation time. 

We use a differential technique to determine the time 
dependence of the pair recombination. After setting the 
pump-probe delay (t) a square wave signal around 100 Hz 
is fed to the laser electronics making the laser pulses hop 
between two time locations. The time separation {5t) be- 
tween these two locations can be set to any value up to 
600 ps. The modulation signal then serves as the refer- 
ence for lock-in detection of the far-infrared signal. The 
signal is the difference in the transmission of two excited 
states occurring at pump to probe delays of t — St/ 2 and 
t + St/2. One can then vary the delay t over the range 
where there is a signal. Subsequent integration gives the 
full relaxation curve at a given temperature. Typical 
results are shown in Fig. ^c). Repeating the measure- 
ment at various temperatures allows us to determine the 
thermal evolution of the signal amplitude and effective 
relaxation time. These two quantities, however, can also 
be determined measuring only the maximum (positive) 
and minimum (negative) points in the derivative signal. 
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FIG. 2: (Color online) Ratio of superconducting-state to normal-state far-infrared transmission for several temperatures. The 
symbols are the experimental data and the lines are strong-coupling BCS fits using the parameters given in Table Q Note that 
the frequency scale for the NbN film is broader than for the other samples. 



Material 


Rn{n) Tc (K) 2Ao( cm"') 


2Ao /ksTc 5S{n 


-^cm-^) 


r^/rS / (nJ / pulse) 




Pb 


70 


7.2 


22.5 


4.5 


0.5 


2.53 


0.5 


0.03 (3.7 K) 


Pbo.75Bio.25 


110 


8.0 


28 


4.9 


0.7 


2.2 


0.5 


0.24 (3.0 K) 


Nb 


61 


5.4 


14 


3.7 





13.3 


2 


0.03 (3.5 K) 


Nbo.5Tio.5N 


110 


9.8 


26.5 


3.9 


0.2 


6 


2 


5 (2.0 K) 


NbN 


70 


13.5 


35.0 


3.7 





1.53 


2 


14 (2.2 K) 



TABLE I: Parameters for Mattis-Bardeen fits of the far infrared superconducting to normal transmission ratio Ts/T„. In all 
cases we used the dirty limit (1/r — > (x>) approximation. The time ratio obtained from the non-equilibrium fits, the laser 
fluences (filling the 5 mm sample aperture) utilized and the excess QP fraction Nqp/Nr are also shown. For the latter the 
number in parenthesis is the lowest sample temperature reached. 



These points are indicated in Fig. ^c) by 'A' and 'B', 
respectively. In a linear approximation, the amplitude of 
the integrated decay curve [squares in Fig.^c)] is propor- 
tional to the maximum slope (point 'A'). This amplitude 
is related to the number of excess QPs {Nqp) through 
curves like the one in Fig.^b), yielding 



Ngp{T) <xAx 



A^(0) Ts(0) 



(1) 



Using the synchrotron pulse width, tyuv, (which deter- 
mines the time resolution) and assuming that the decay 
is represented by a single exponential (at least in the be- 



ginning of the relaxation), the effective relaxation time is 
given by 

A 

^eff = -tvUV X 



(2) 



The photoinduced spectra were also determined in a 
differential mode. As the lamellar spectrometer works 
in a step-scan mode, the dithering method described in 
the previous paragraph can be used straightforwardly. To 
take advantage of the fast scan mode in the Bruker inter- 
ferometer we could not utilize a lock-in based detection 
system. To obtain the differential spectra we recorded 
a set of 33 spectra alternating the two laser delays {t 
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and t + 5t). Each spectrum has n scans except for the 
first and the last spectra which are taken at the same 
delay with n/2 scans. All the measurements at the de- 
lay t are summed up in T and all those taken at t + St 
are added into T'. We can then determine the differen- 
tial spectrum as — AT/T = (T — T') /I. This procedure 
minimizes the influence of the synchrotron light intensity 
variations arising from storage ring current decay. 

As final remark, we are interested in the low-fluence 
regime. To verify that we are in this limit we set the 
pump-to-probe delay so that the differential signal is 
maximum [point 'A' in Fig. ^c)]. A necessary but not 
sufficient condition for a system to be in the low-fluence 
regime, is to have this signal varying linearly with the 
number of incident photons. By changing the laser power 
we can determine the boundary of the linear regime. This 
procedure was done for each sample at the lowest temper- 
ature measured, where the number of thermally produced 
QPs is the lowest. 

III. THE EQUILIBRIUM FAR-INFRARED 
RESPONSE 

The transmission of light thorough a thin film and into 
a substrate of refractive index n is given in the long- 
wavelength limit by 



y = yi + iy2 — Zoad is the film's dimensionless complex 
admittance with a being the optical conductivity, d the 
film thickness and Zq « 377 O the vacuum impedancei^ 
Equation 13 is valid when the film thickness is much less 
than the penetration depth and the measurement wave- 
lenght. 

The normal-state frequency-dependent infrared re- 
sponse of metals and alloys is fairly well described by 
the Drude model. The electrodynamics of BCS super- 
conductors was calculated by Mattis and Bardeen^ for 
dirty limit superconductors (l/r 3> 2 A) and an exten- 
sion of those calculations was proposed by Zimmermann 
et almL for an arbitrary scattering rate. All these ex- 
pressions, valid in the BCS weak coupling limit, give 
the ratio (trf -I- icrf ) /ctq, where is the superconduct- 
ing state complex conductivity and ao is the normal- 
state dc conductivity. The theory for the strong-coupling 
superconducting-state optical conductivity was produced 
by Nam4^ Ginsberg et al.'^^ pointed out that the low fre- 
quency approximation for these expressions corresponds 
to a decrease in the condensate weight. This decrease 
manifests itself as the subtraction of a factor 2 SS/ttuj 
from the imaginary part of the Mattis-Bardeen conduc- 
tivity. 

Note that Eq. 13 depends only on the product ad and 
not on the individual values of these two quantities. 
From the normal state transmission one can determine 
straightforwardly the sheet resistance Ru = 1/ u^d which 



can, in turn, be used to calculate the absolute values of 
d. This procedure works rather well for our samples 
because they are in the dirty limit, making lj <^1/t over 
the entire far infrared, so that the normal-state optical 
conductivity is essentially identical to the dc conductiv- 
ity. 

We show in Fig. 13 the ratio Ts/Tn for all our sam- 
ples. (Subscripts or superscripts S and N correspond 
respectively to the superconducting and normal states.) 
The parameters used in the fits are given in Tabled In 
each sample we fitted the lowest temperature curve by 
varying the gap and the strong coupling coefficient 5S. 
The experimentally determined temperature was then 
used to obtain the Mattis-Bardeen behavior (without 
changing any other parameter) for the remaining sam- 
ple data. Within experimental errors, the fitted trans- 
missions agree with the data. The fitting parameters 
are also in agreement with known values for these ma- 
terials. Hence, the Mattis-Bardeen formalism and also 
our approximation in Eq. El describe correctly the super- 
conducting response of our films justifying their use in 
the analysis of the photoinduced state fSec lfVC)) and in 
evaluating the sample temperature (see the Appendix). 



IV. THE EXCESS QUASIPARTICLE STATE 
A. Theory 

Figure|21(a) shows the density of states at K of a BCS 
superconductor. Light having a photon energy hv > 2Ao 
is absorbed and breaks a Cooper pair, producing two QPs 
above the Fermi energy. In our experiments the pho- 
ton energy is much higher than the superconducting gap 
and thus the QPs are created far above the Fermi level 
with an energy ~ hv /2. Very quickly (in sub-picosecond 
times) these high-energy QPs relax via electron-electron 
and, eventually, electron-phonon scattering. These re- 
laxation processes break more pairs, and the QPs quickly 
settle to energies close to A. In such a cascading process, 
each QP initially created by the laser pulse will generate 
^ hv/2/S. QPs close to the Fermi level. 

A similar branching process occurs in the phonon sec- 
tor. Phonons of energy higher than 2A can break pairs, 
relaxing to 2A whereas phonons of energy less than 2A 
do not participate in the pair breaking process. There- 
fore, the system eventually reaches the state depicted in 
Fig. Elb). This state is characterized by an excess (with 
respect to thermodynamic equilibrium) of unpaired QPs 
of energy A and, as shown by Owen and Scalapino,"'^^ by 
a reduced energy gap. 

The effective relaxation process must take into account 
a phonon bottleneck effect, shown schematically in Fig. 
|3Jc). Two QPs having energy A take a time Tfi to re- 
combine into a Cooper pair. At recombination, a 2A 
phonon is emitted to carry the excitation energy of the 
QPs. The phonons with energy 2A can break a Cooper 
pair, in characteristic time tb, creating two QPs at the 
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FIG. 4: (Color online) (a) Effective relaxation time (tc^ ) cal- 
culated from the intrinsic QP recombination time (ta), the 
pair breaking time by phonons (tb) and the phonon escape 
time (r^) using Eqs. [HI |H1 and ||l with = r]^ = 100 ps. 
r-y = 300 ps is taken to be temperature independent, (b) 
Temperature dependence of the excess number of QPs from 
Eq. 111! The dashed line was calculated using the same pa- 
rameters as those used in panel (a). For the other curves, the 
number next to each line indicates the ratio t^/t£ utilized. 



FIG. 3: (Color online) Ratio of the superconducting (Ns) and 
normal (A''jv) densities of states at K for a BCS supercon- 
ductor (a) in equilibrium and (b) after photoexcitation and 
initial cascading effects but before QP recombination. The 
broken pairs in the photoexcited state produce the empty 
states below Ep and occupied states above. Consequently, the 
superconducting gap is depleted. Panel (c) illustrates the re- 
laxation process with competition between QPs and phonons. 



gap edge. A steady-state dynamic balance is established 
between the phonons and QPs. Eventually, however, the 
2A phonons become depleted; either they relax anhar- 
monically to lower energies or they leave the film. (The 
phonons can escape into the substrate or into their sur- 
rounding environment, such as the helium bath.) The 
phonon escape time is denoted t^. 

The general case for any density of excess QPs has been 
discussed by Rothwarf and Taylor.^ We are interested in 
the low-fluence limit, where the number of QPs created 
by the laser light is much smaller than the thermally 
broken pairs. In this case the differential equations can 
be linearized givingiAft 



IP 



2N, 



qp 



dt 



TB. 



2Nn 

TB 



(4) 



and 



dt 



TB 



N, 



1111 + 1111 ] + '_JiP, (5) 

- TR 



Nqp is the number of excess QPs at energy A and Nn 
is the number of phonons with energy 2 A. Assuming 
exponential decays for Nqp and Nn, one can solve the 
system given by Eqs. 0] and |S1 The effective relaxation 
(Teff) time is given by: 



1 

Teff 



2t; 



1- Wi- 



(27 



-1^2 



(6) 



Another solution with the positive sign before the 
square root gives the time for phonons and QPs pop- 
ulations to reach equilibrium. Figure 0] (a) shows that 
for temperatures above ^ OATc but not very close to Tc 
the phonon escape time is much larger than tr or tb ■ 
Hence Eq. can be approximated by: 



TR 

2tb 



(7) 



In our numerical analysis we utilized the full expression 
given by Eq. However, the temperature range where 
Eq. is valid corresponds to the region covered by most 
of our experiments and we will use this expression for 
some qualitative analysis. 

The phonon escape time depends on several factors 
such as the sample quality, the film thickness, acoustic 
mismatch between film and substrate, and the presence 
or absence of liquid helium surrounding the film.^'* There 
is no detailed calculation for the temperature dependence 
of Tj. However, at temperatures much smaller than the 
Debye temperature, the Kapitza resistance in a material 
interface behaves as T^.'^^ As a higher Kaptiza resistance 
makes it harder for phonons to leave the material, in a 
first approximation we can assume that =T° + aT^. 
Several scattering rates for metallic superconductors were 
calculated by Kaplan et al^ In particular, tr and tb can 
be obtained from intrinsic materials parameters: 
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FIG. 5: (Color online) Relative variation in the overall far-infrared transmission — AT/T as a function of the delay between 
pump and probe pulses. 



and 



TRiiu) 1 + fiuj) \kBTj 



+s ^{n-iu)^-d^ 



[nin) + l]fin-u;) (8) 



1 1 



[1 - fin) - fico - n)] 



(9) 



Here, t° = hZi{0)/2Trb{hTcf and t° 
fi7V/47r^iV(0)(a^) Aq are respectively the intrinsic 
pair recombination time close to Tc/2 and the phonon 
pair-breaking time at K. All energies in Eqs.|HlandEl 
are measured in units of Aq. S = A/Aq is the reduced 
gap, Zi{Q) is the QP renormalization factor, A^(0) is the 
electron single particle DOS at Ep, N is the density of 
ions, and (a)^ is the electron-phonon coupling function 
averaged over the whole phonon spectrum. In Eq. |Slthc 
Eliashberg function a'^F{n) is approximated by bil^, 
but the correct function can be used if it is known, lu 
is the energy of QPs just before recombination or the 
energy of photons available for pair breaking. n{fl) and 
/(fi) are the Bose and Fermi factors, respectively. 

Figure EJa) shows tr and tb calculated from Eqs. |H1 



and El assuming a weak coupling BCS temperature de- 
pendence for the gap, QPs with energy A and phonons 
with energy 2.1 A. We note that diverges at K while 
Tb remains finite. Therefore, at low temperatures, all 
the photon energy absorbed in the pair breaking process 
will appear as excess QPs. At higher temperatures the 
photon energies will end up distributed between QPs and 
phonons so that 

Ngp{0)Ao = Ngp{T)AiT) + NniT)2A{T). (10) 

When the phonon and QPs populations are in equilib- 
rium, Nqp/rii w Nq/tb and, thus, 

N,p{T) _ Aq 1 

NgpiO) A{T)1 + 2tb/tr ^ ' 
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With Eqs. (or|7I) and II II we can describe both the ex- 
cess number of QPs and the effective relaxation time as a 
function of the temperature. Simulations of the effective 
relaxation time and number of excess QPs are shown in 
Fig. 2] Note that whereas the effective relaxation time 
depends on the phonon escape time (a parameter that 
varies from sample to sample), the excess number of QPs 
depends only on the ratio t^/t''^. In addition, its strong 
dependence on t^/t^ allows for an accurate determina- 
tion of this ratio. 



B. Integrated response 

Figure shows the overall photoinduced far-infrared 
transmission for five films as a function of the time be- 
tween pair breaking laser pulse and far-infrared probe 
pulse. For the samples where the data have been taken 
on a cold finger cryostat (Pb and Pbg. 75610.25) there is an 
initial fast relaxation followed by a long-lived tail. Data 
taken with the samples in contact with liquid or gaseous 
He [Nb, Nbo.5Tio.5N and NbN] do not show this long- 
lived tail. However, this tail reappears in the data for the 
Nb sample when measured on the cold- finger cryostat. 
We therefore assign this long relaxation time process to 
film and substrate thermal effects rather than to intrinsic 
properties of the superconducting films. Heating issues 
are very important in pump-probe measurements and are 
discussed thoroughly for our samples in the Appendix. 

Disregarding the long-lived tail in the Pb and 
Pbo.75Bio.25 samples, the relaxation back to equilibrium 
can be described in all samples and at any temperature 
by a single exponential. We can then use the approxi- 
mations given by Eqs.^andl^to study the temperature 
dependence of the excess number of QPs and the effective 
relaxation rate. 

Figure shows the effective relaxation time as a func- 
tion of temperature for the Nb and Nbo.5Tio.5N samples. 
The fits describe the data properly at temperatures above 
Tc/2 but fail below this temperature, with the measured 
times being smaller that the calculations. One possible 
explanation for the discrepancy relates to the QP energy 
before recombination. In our calculations we assumed 
that the QPs were sitting at an energy A prior to recom- 
bining into pairs and we disregarded recombination of 
higher energy QPs. Calculations using Eq. Oshow that 
tr decreases for increasing w. Indeed, it is possible to 
improve the fits between 0.2 and 0.6 Tc in Fig. |H|by in- 
creasing uj. This happens, however, at the expense of the 
fitting quality at higher temperatures. A complete de- 
scription of this effective time may require recombination 
of QPs sitting at a temperature dependent distribution of 
energies. Another possibility is the effect of the smaller 
number of thermally broken pairs at low temperatures 
pushing the system into the high-fluence regime: the in- 
creased number of unpaired QPs available decreases the 
recombination time. 

The data in Fig.[7|was calculated using Eq.^and then 
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FIG. 6: (Color online) Effective relaxation time for Nb and 
Nbo.5Tio.5N. The points are the experimental data trans- 
formed using Eq. |21 and the solid lines are simulations using 
Eq.|S| In both cases a small contribution to is necessary 
to describe the rising part of the data above 0.8 T /Tc. The 
error bars are obtained from statistical standard deviations 
from approximatelly 500 measurements at each temperature. 




FIG. 7: (Color online) Thermal dependence of the number of 
excess QPs in all our samples. The solid lines are fits using 
Eq.ini 



corrected using curves similar to the one in Fig. Q^b). 
The fits for Tf.ff{T) and Nqp{T) do not allow for a de- 
termination of the absolute values of andr^ as Eqs.[3 
and depend only on the ratio of these two time con- 
stants. This ratio however, is representative of the strong 
coupling character of the system. As shown in Fig. |H1 a 
long phonon pair breaking time tb compared to cor- 
responds to a large 2Ao/A:bTc ratio. 
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FIG. 8: (Color online) lAo/ksT^ plotted against tr/tb- The 
solid symbols are the results of this work. The open symbols 
are obtained from Kaplan et al^ calculations. The dashed 
line is a guide for the eye. 



C. Photoinduced spectra 




1 - 
- 

-1 '- 



-1 - 



I I I 



I I I 




(Nb,Ti)N 



I I I 



28A = -0.1 cm 
I I I I I I I I I I I I I I I 
PbBi 
3.0 K 

28A = -0.12cm'' - 
I I I I I I I I I I I I I I I 
Nb 
3.5 K 

25A = -0.23 cm'' J 
I I I I I I I I I I I I I I I J 
2K 25A = -0.26 cm ' 

8 K 2SA = -0.05 cm"' 
I I I I I I I I I I I I I I I ~ 




NbN - 2.2 K 




28A = -0.1 cm J 
I ■ ■ ■ I ■ ■ ■ 



20 



40 



60 



80 



100 



Frequency (cm" ) 



The measurements discussed up to now were all inte- 
grated photoinduced transmission; no spectral features 
can be studied in this approach. Although those mea- 
surements determine the relation between the character- 
istic system relaxation times, one obviously cannot get di- 
rect information on the gap depletion predicted by Owen 
and Scalapino.^'* By spectroscopically resolving the probe 
pulses we can compare the excited state and equilibrium 
spectra and directly measure the superconducting gap 
weakening. 

In this case, we are mostly interested in the photoin- 
duced response with pump and probe in coincidence. Al- 
though one can take measurements at several pump-to- 
probe delays, the time dependence of the gap depletion 
follows the effective relaxation time T^g, already deter- 
mined in the integrated measurements. Using the present 
synchronized method to measure the photoinduced spec- 
tra has the advantage of minimizing long-lived lattice 
heating effects and detecting the true changes produced 
by the excess QPs state. 

Equation O shows that the superconducting-state 
transmission depends on the optical conductivity. A 
slight changes in a corresponds to changes in the trans- 
mission as: 



AT _ 7{c7 + Sa) - 7{a) 



(12) 



where da = a{A + 6 A) — (t(A). SA is the change in 
the gap arising from photoexcitation. In this picture, 
the equilibrium response Mattis-Bardeen equations can 
be used to describe a photoinduced state characterized 
by a gap weakening. 

Figure 1^1 depicts a few AT/T spectra around the gap 
energy using the pump fluences shown in Table ^ The 



FIG. 9: (Color online) Photoinduced spectra (pump and 
probe in coincidence) for all materials studied. The symbols 
are the experimental data and the lines are fits using Eg. 1121 
and Mattis-Bardeen optical conductivities. 



solid lines are simulations using Eq. ^] but keeping the 
same parameters used in the fits for Ts /7n- These differ- 
ential spectra allow us to determine directly the gap shift 
(2 5A) produced by the excess QPs state. The magni- 
tude of 2 SA controls the height of the jump around 2A 
where —AT/T changes sign. The BCS fits are able to 
describe the data above 2A and the size of the —AT/T 
jump. Conversely, the quality of the low-energy fit is gen- 
erally poor, with the data not showing as sharp a discon- 
tinuity as predicted by the theory. However, the basic 
physics appears to be what is expected: the supercon- 
ducting gap is shifted to lower energies in the presence 
of the laser light, as the superfluid density is depleted 
and the zero-frequency superconducting delta function is 
weakened. 

The knowledge of the gap shift also allows us to verify 
whether or not our data were in the low-fluence regime. 
In BCS theory, the superfluid concentration Ug is of the 
order of A for T < and A^ for T « T^. The ratio SA 
then tells us how many excess QPs we have in the system 
and this number can be compared to 1 — A^(0)/A^(T) 
representative of the number of thermally broken pairs 
Nt- The low-fluence limit is then given by the condition 



Nt 



5A 
~A 



A2(T) 



< 1. 



(13) 



Values for Ngp/Nx are shown in Table HI at the low- 
est temperature measured in each sample. We no- 
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tice that the low-fluence hmit is indeed fulfilled for Pb, 
Pbo.TsBig 25 and Nb at all temperatures but not for low 
temperatures in Nbo.5Tio.5N and NbN. This effect may 
then explain the failure of our simple description of t^jj at 
low temperatures and the better fits for — AT/T obtained 
in Pb, Pbo.75Bio.25 and Nb. It is also compatible with a 
poor — AT/T fit at low temperatures for Nbo.5Tio.5N and 
a good fit for the same sample at 8 K. Indeed, at this 
temperature we calculated that Nqp/Nr ~ 10~^. The 
understanding of the high fluence regime requires more 
data but deviations from our observations here are ex- 
pected. In cuprate superconductors, for example, it was 
shown that the recombination dynamics is strongly de- 
pendent on the number of broken pairs ^ 



V. CONCLUSIONS 

We measured the equilibrium and photoinduced 
time dependent far-infrared response of several Pb, 
Pbo.75Bio.25, Nb, Nbo.5Tio.5N and NbN films in the su- 
perconducting state. We determined the effective relax- 
ation time and number of excess quasiparticles created 
after the system has pairs broken by a fast laser pulse. 
Typical relations in the nanosecond range were observed. 
The most relevant quantity to characterize the relaxation 
process is the ratio between the intrinsic pair recombina- 
tion time and the pair breaking time by phonons. The 
experimental results are in reasonable agreement with 
the linearized Rothwarf and Taylor^ relaxation equations 
proposed by Gray.^*' The gap shrinkage expected for the 
excess quasiparticle state was directly observed in the 
photoinduced far-infrared spectra. From a detailed anal- 
ysis of the gap shift we can monitor the regions where 
the low-fiuence limit approximation is valid. 
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APPENDIX: HEATING EFFECTS 

One of the major experimental issues in all optical 
pump-probe spectroscopy is the ability to separate the 
intrinsic photoinduced effects from those produced by 
lattice heating. In our measurements, we observed two 
distinct lattice heating effects. The first is an average 




Frequency (cm"') Delay (ns) 

FIG. 10: (Color online) (a) Ts/Tjv for a Pb film as a function 
of laser power in CW mode measured at a nominal tempera- 
ture of 4.77 K. The solid lines are BCS fits assuming higher 
eff'ective lattice temperatures due to laser absorption, (b) 
Efi'ective relaxation times for 3 different Pb films. The fast 
initial relaxation is the same for all films (the slight difference 
is due to different synchrotron pulse widths setting the time 
resolution). The long lived tail changes randomly from film to 
film and we assign it to different thermal contact conditions. 



'static' heating produced by long term photon absorp- 
tion. The second is a quick lattice heating that appears 
when the system is not able to dissipate the photon en- 
ergy fast enough. 

To evaluate and almost completely eliminate the 
'static' effect, we can use the changes in the far-infrared 
spectrum. This is done by setting the laser to its CW 
mode and measuring the far-infrared transmission as a 
function of laser power. In this case the laser pulses are 
at random times and therefore not synchronized to the 
infrared pulses. Alternatively one can also set the delay 
between pump and probe to a time much larger than the 
system excess QP relaxation time. In our samples both 
methods yielded the same results. Figure [rUT a) shows a 
series of Ts/Tat spectra for a Pb sample taken at iden- 
tical temperatures but using different laser powers. The 
increasing laser power changes the far-infrared transmis- 
sion which can be described by a Mattis-Bardeen con- 
ductivity with a higher effective temperature. By fitting 
these spectra we can infer the temperature change pro- 
duced by the laser and search for the laser power that 
minimizes this change. In our measurements, we made 
sure that this 'average' lattice heating is less than 0.5 K. 
It is very important to remark that this effective lat- 
tice temperature is not the same T* discussed by Parker 
and Williamsi^. In their formalism T* is an effective 
electronic temperature produced by an excess number of 
QPs. In our heating determination, the origin of the ef- 
fective temperature is at the heat produced by the decay 
of electron hole-pairs created by the laser in the film. 

The second heating effect is harder to discriminate 
from intrinsic processes. While discussing the decays for 
Pb and Pbo.75Bio.25 in Fig.|31we assigned a long lived tail 
on the time dependent — AT/T to lattice heating. We 
have two major motivations to make this statement. As 
these films have been grown in situ at 80 K and needed 
to be kept cold at all times, we were constrained to mea- 
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sure them in a cold finger cryostat. The heat dissipation 
in this case only passes through the area in contact with 
the cryostat tip. As we are measuring the sample trans- 
mission, this thermal contact does not happen with the 
area under illumination. As shown in Fig. llOf b). differ- 
ent Pb samples show the same short decay but the long 
lived tail changes from sample to sample indicating an 
extrinsic process. Because the sample thermalization de- 
pends on the thermal contact between film and substrate 
and between substrate and cold finger, and because these 
factors were not always identical, this long lived tail is 
likely to be the signature of the lattice heating. The 
second motivation to make this assignment comes from 
measurements on a Nb film. The relaxation curves for 
this sample was measured in both cryostats; the long tail 



was only present for measurements where the sample was 
in vacuum in a cold- finger cryostat. When the sample is 
in He gas or liquid, the area under illumination is being 
cooled down by the helium, improving the heat dissipa- 
tion. 

In summary, we can easily eliminate the 'average' 
static sample heating by carefully checking the appropri- 
ate laser power used as a pump. A second, time depen- 
dent heating effect, can be identified by comparing the 
decays in different samples. This time dependent heating 
effect has not been observed in measurements made with 
the sample in contact with He liquid or gas and indicate 
that it originates in a bad energy dissipation through the 
cryostat cold finger. 
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